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Abstract: The Mad2 protein, with two distinct conformations of open- and closed-states, is 
a key player in the spindle checkpoint. The closed Mad2 state is more active than the open 
one. We carried out conventional and targeted molecular dynamics simulations for the two 
stable Mad2 states and their conformational transition to address the dynamical transition 
mechanism from the open to the closed state. The intermediate structure in the transition 
process shows exposure of the p6 strand and an increase of space around the binding sites 
of P6 strand due to the unfolding of the P7/8 sheet and movement of the P6/4/5 sheet close 
to the aC helix. Therefore, Mad2 binding to the Cdc20 protein in the spindle checkpoint is 
made possible. The interconversion between these two states might facilitate the functional 
activity of the Mad2 protein. Motion correlation analysis revealed the allosteric network 
between the pi strand and P7/8 sheet via communication of the P5-aC loop and the P6/4/5 
sheet in this transition process. 

Keywords: molecular dynamics simulations; open/closed Mad2; dynamical transition 
mechanism 
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1. Introduction 

The spindle assembly checkpoint (SAC) proteins monitor a series of events in mitosis to ensure the 
accuracy of chromosome segregation prior to cell division, and are indispensable for chromosome 
stability [1-4]. Mitotic arrest deficiency 2 (Mad2), which is a critical molecular component of the SAC 
proteins, is able to inhibit the anaphase-promoting complex or cyclosome (APC/C), combining with 
Cdc20 (an activator of APC/C) to delay separation of premature sister-chromatid through the binding 
of Mad2 to Cdc20 [5-14]. Premature sister-chromatid separation leads to aneuploidy (abnormal 
number of chromosomes), which contributes to cancer progression [5,15-18]. Mad2 binding to Cdc20 
is helpful for accurate chromosome segregation to prevent cancer development. Therefore, it is 
important that the active Mad2 protein binds to Cdc20 in the mitotic process. 

Usually, a given protein has a unique and native three-dimensional structure with the lowest 
possible Gibbs free energy [19]. However, Mad2 shows the characteristics of a "metamorphic" protein, 
and is able to adopt two distinctly different conformations, i.e., the open and closed Mad2 states under 
native conditions [20-22]. The first high-resolution X-ray crystal structure of the open Mad2 state was 
determined by Luo et al. in 2000 [23]. Since then, significant progress has been made in filling the 
structural space of this family. Luo et al. again reported the structures of several closed Mad2 states in 
2002 and 2004 [20,24]. Subsequently, other crystal structures of the Mad2 monomer and dimer were 
also reported [14,25,26]. These crystal structures revealed that the structure of the closed Mad2 state is 
strikingly different from the open Mad2 state. In the open Mad2 state, the structure of Mad2 consists of 
three layers: a central layer formed by three a-helices (aA, oB and aC), a large six P-sheet ordered as 
P7-P8-P6-P4— p5-pi in the fi-ont layer, and a short P-hairpin (P2, P3) in the back layer. In the closed 
Mad2 state, the predominant differences in features of this state with respect to the open one are the 
large five P-sheet ordered as P6-P4— P5-P8"-P8' and a A^-terminal pi strand forming an additional helix 
attached at the aA helix. Moreover, the P6 strand is exposed in the closed Mad2 state, while it is 
blocked in the open Mad2 state. It has been reported that Cdc20 easily forms an edge-on interaction 
with the P6 strand at the ligand-binding site in the closed Mad2 state. The closed Mad2 monomer or 
dimer is more functionally active than the open one, which might facilitate Cdc20 binding in order to 
impart fidelity chromosome segregation. 

Cytosolic Mad2 in human cells is mainly represented by an open Mad2 monomer [20]. Then, the 
active closed Mad2 conformation should be initially transformed by the open Mad2 monomer with or 
without the Madl-Mad2 template [27-31]. Therefore, the conversion of the open state to the closed 
state in the Mad2 protein plays a key role for Cdc20 binding and has become an area of considerable 
interest over the past decades. It is well-known that the open Mad2 state is able to spontaneously 
convert to the closed Mad2 state [20]. Based on the difference between the open and closed states, 
their conversion should involve the transformation of the P7/8 sheet and pi strand. Some experimental 
data suggests that the open Mad2 state undergoes a conformational change to the closed state via a 
high-energy intermediate conformation, which might favor Cdc20 binding [26,32]. On the other hand, 
the theoretical investigations on the Mad2 protein are very limited so far. Therefore, it remains unclear 
what the dynamical transition mechanism of the two states is and what the exact allosteric network in 
this conversion is. Because it is difficult to experimentally isolate the intermediate structures during the 
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conformational transition process, it is highly valuable to explore the conformational transition process 
in more detail by performing computational simulations at the atomic level. 

In the present study, we carried out conventional molecular dynamics (CMD) and targeted 
molecular dynamics (TMD) simulations in order to explore the conformational transition pathway in more 
detail. The results obtained from the simulations provide valuable insights into this conformational 
transition process with respect to how one state gradually transforms into another. 

2. Results 

The initial structures of the open and closed Mad2 protein states (denoted as 0-Mad2 and C-Mad2, 
respectively, for convenience) used in the CMD simulations were based on previously reported and 
deposited X-ray crystal structures (PDB codes IDUJ and 2VFX (K chain)) [14,23]. The 0-Mad2 
model consists of three a helices (aA: 18-37, aB: 61-75, and aC: 122-142) and eight P sfrands (pi: 
11-13, P2: 43^9, P3: 53-58, P4: 81-90, p5: 97-105, P6: 150-160, P7: 170-174 and P8: 184^191). 
The missing residues (Metl-GlylO, Metl96-Asp205) in this model were built using the loop 
search method in the Swiss-Pdb Viewer (http://spdbv.vital-it.ch/). The C-Mad2 model consists of 
three a helices (aA: 12-37, aB: 61-75, and aC: 122-142) and seven p sfrands (P2: 43^9, P3: 53-58, 
P4: 81-90, P5: 97-105, P6: 150-160, P8': 177-188, and P8": 190-200). The mutated residue Alal3 in 
the C-Mad2 model was modeled by starting from the original residue Leul3. 

To explore the detailed conformational transition pathway, we first carried out CMD simulations on 
the 0-Mad2 and C-Mad2 states to understand the stability of the two Mad2 structures, followed by 
TMD simulations associated with the transition from the open state to the closed state and vice versa to 
explore the dynamic transition between the two states of the protein. Depicted in Figure 1 are the plots 
for the time-dependent root mean square deviations (RMSDs) for the positions of all backbone atoms 
in the CMD-simulated 0-Mad2 and C-Mad2 structures, revealing that the CMD trajectories tended to 
be stabilized after about 5 ns; the corresponding energies also tended to be stabilized after about 5 ns. 
Therefore, for each protein state, trajectory analysis was performed to extract the equilibrated 
conformation between 30 and 50 ns simulation time, recording 1 0,000 snapshots (one for each 2 ps) 
for each frajectory. Depicted in Figures 2 and SI are the plots for the RMSDs of the positions of all 
backbone atoms in the TMD-simulated protein structures from the corresponding final structures 
associated with the two TMD frajectories for the 0-Mad2<-^C-Mad2 conversion. The plots depicted in 
Figures 2 and S 1 reveal that the backbone atoms in the two TMD simulations all reached the target 
conformations within 12 ns. The detailed structural fransition results obtained from the CMD and 
TMD simulations are discussed below in detail. 
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Figure 1. Root mean square deviation (RMSD) values of all backbone atoms of the 
0-Mad2 and C-Mad2 models with respect to the corresponding starting structure for the 
simulations of 0-Mad2 and C-Mad2 models. 
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Figure 2, Root mean square deviation (RMSD) values of all backbone atoms and all atoms 
of the whole protein are shown for the 0-Mad2^C-Mad2 transition. 
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2.1. Stability of the 0-Mad2 and C-Mad2 States 

Based on the CMD simulations (see Figure 1), we obtained two dynamically stable states (the open 
and closed states) of Mad2. For each state, the average structure extracted from the CMD trajectory is 
shown in Figure 3. According to our analysis on the simulated structures, both the 0-Mad2 and 
C-Mad2 structures simulated reflect the characteristics of the X-ray crystal structures very well. 
For example, these two states consist of three layers: a central layer formed by three a-helices, a large 
six- or five-stranded P-sheet in the fi-ont layer, and a short P-hairpin in the back layer. Specifically, in 
the 0-Mad2 state, the three a-helices in the central layer are aB-aA-aC, a large P7-p8-p6-p4-p5-pi 
P-sheet in the Iront layer, and a short P2-P3 P-hairpin in the back layer arranged orderly fi-om left to right 
in Figure 3a. In the C-Mad2 state, except for the large P6-P4-P5-P8"-P8' P-sheet rearranged in the front 
layer, similar components in the central and back layers can be found in Figure 3b. The most 
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remarkable differences between the 0-Mad2 and C-Mad2 states are: (1) the P7/8 sheet is located on 
the left side of the p-sheet in 0-Mad2 and the p878' sheet located to the right in C-Mad2; 
(2) the pi strand in 0-Mad2 forms a new helix attached to the A^-terminus of aA in C-Mad2. 

Figure 3, MD-simulated Mad2 structures for the open state (a) and closed state (b). 



2.2. The Conformational Transition from 0-Mad2 to C-Mad2 

Based on the structural features of the metamorphic protein Mad2, we performed TMD simulations 
to explore the conformational transition process between 0-Mad2 and C-Mad2. Based on the two 
TMD simulations for the 0-Mad2-<-»^C-Mad2 conversion, the forward and reverse transitions almost 
followed similar structural conversion processes computationally (see Figures 2, 4, SI and S2). 
Consequently, we only present the conformational transition from 0-Mad2 to C-Mad2 starting from 
the equilibrated 0-Mad2 structure. During the transition process, three transition conformations 
(I (-3.5 ns), II (-6.5 ns). III (-9.5 ns)) of the Mad2 protein were extracted from the TMD simulation 
trajectory and were simulated via the CMD method as initial structures. The corresponding 
time-average structures are shown in Figure 4. In the first step of the transition (from 0-Mad2 to I), the 
pi strand in the 0-Mad2 state extends to a line loop and, simultaneously, the PV/S sheet unfolds to 
form conformation I with the destruction of their hydrogen bonds and hydrophobic interactions. In the 
second step of the transition (from I to II), the unfolded pi strand in conformation I traverses through 
the P5-aC loop close to the aA helix and forms a new helix attached to the A^-terminus of aA helix to 
produce conformation II, which causes the large fluctuation of the P5-aC loop and favors the P7/8 shift 
through emptying the right region of the P6/4/5 sheet. Simultaneously, the P6/4/5 sheet in the front 
layer moves close to the aC helix in the central layer (see Figure 4), which results in the exposure of 
binding sites around the P6 strand and an increase in the space between the P6 strand and the unfolded 
pV/S sheet, ultimately increasing its activity to combine with the Cdc20 protein. Because the 
structural characteristics of II showed that the new helix, attached to the A^-terminus of aA helix is 
completely formed and that the P6/4/5 sheet is closest to the aC helix in the central layer 
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(see Discussion Section), II was regarded as the intermediate during the conversion of 0-Mad2 to 
C-Mad2. In the third step, i.e., the conformational transition from II to III, the unfolded P7/8 sheet in 
conformation II moves from the left side of p6/4/5 sheet to its front to produce conformation III. In the 
fourth step of the transition (from III to C-Mad2), the unfolded P7/8 sheet in conformation III moves to 
the right side of the P6/4/5 sheet to form a new p8"/8' sheet and produces the C-Mad2 state through 
formation of new hydrogen bonds and hydrophobic interactions (see Figure 4). In summary, the 
transition process from 0-Mad2 to C-Mad2 involves the initial unfolding of the P7/8 sheet followed by 
the transfer of the unfolded P7/8 sheet. Moreover, the pi strand traversing through the P5-aC loop is a 
prerequisite for P7/8 translocation. The intermediate II is characterized by the exposure of the 
P6 strand and the increase in space of binding sites around P6, which is possibly a favorable 
conformation to combine with the Cdc20 protein to mediate APC/C inhibition. 

Figure 4. The three-dimensional structures of the conformational transition pathway 
determined by TMD from 0-Mad2 to C-Mad2. 



2.3. Correlation Analysis along the Transition Pathway from 0-Mad2 to C-Mad2 

To explore the allosteric communication of position changes for the pi strand and the p7/8 sheet 
along the transition pathway from the 0-Mad2 state to the C-Mad2 state, motion correlations of 
Ca atoms were analyzed and the corresponding cross-correlation map constructed from the trajectory 
is given in Figure 5 for the conversion of 0-Mad2^C-Mad2. The results showed that the motion 
correlations between the residues ranged from high anticorrelation (blue) to high correlations (red). 
The results in Figure 5 demonstrate that the large correlated motions of the pi strand vs. the pS-aC 
loop occur remarkably with the large correlated motions of the pS-aC loop vs. the P6/4/5 sheet, which 
predicts that the allosteric communication between the pi strand and the P6/4/5 sheet occurs via the 
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P5-aC loop region. Furthermore, the motions of the P6/4/5 sheet significantly correlate with the 
motions of the P7/8 sheet, which predicts the allosteric communication between the P6/4/5 and the 
P7/8 sheets. 

Figure 5. Dynamical cross-correlation map of the transition from 0-Mad2 to C-Mad2 with 
specific sub-region squares in black. 
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3. Discussion 

3.1. Conformational Characteristics along the Transition Pathway from 0-Mad2 to C-Mad2 

To explore the details of the structural changes along the transition pathway from 0-Mad2 to 
C-Mad2, we used the aA and aC helices as reference points to calculate the variations in mass center 
distances between the pi strand and the aA helix, between the P7/8 sheet and the aC helix, and 
between the P6/4/5 sheet and the aC helix. The corresponding results are shown in Figure 6. It can be 
seen from the conversion of 0-Mad2 to C-Mad2 that the distance between the pi strand and the 
aA helix increases from 11.4 to 22.5 A, which corresponds to the pi strand traversing through the 
P5-aC loop and the emptying the right region of the P6/4/5 sheet, consequently favoring the PV/S 
translocation. Moreover, the distance between the P7/8 sheet and the aC helix decreases from 26.2 A 
in the 0-Mad2 state to 14.8 A in the C-Mad2 state, which represents the movement of the P7/8 sheet 
from the left side of the P6/4/5 sheet to its right side. Interestingly, this distance variation is slower 
than that between the pi strand and the aA helix (see Figure 6), which suggests that the pi strand 
traversing through the P5-aC loop occurs prior to the P7/8 translocation, as discussed in Section 2.2. 
In particular, in the conformational transition from 0-Mad2 to II, the distance between the P6/4/5 sheet 
and the aC helix decreases from 18.1 A in the 0-Mad2 state to 15.0 A in the intermediate II, which 
corresponds to the movement of the P6/4/5 sheet in the Iront layer closest to the aC helix in the central 
layer to form the intermediate. As expected, the moving distance of the p6/4/5 sheet with 3.1 A 
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predicts that the Cdc20 binding sites, i.e., residues Leu61, Tyr64, Asn67, Trp75, Aspl52-Tyrl56, 
Aspl58, Aspl60-Leul61, VaI163, Lysl66-GIul69, Prol72, Phel74 around the p6 strand [24], in the 
intermediate could be considerably exposed compared with the 0-Mad2 state (see Figure 7a). 
To further address the exposed space of Cdc20 binding sites in the intermediate, the CASTp program 
was used to calculate the size and shape of its pocket compared with that in the 0-Mad2 state. Such 
analysis was shown to be able to provide a comprehensive and detailed quantitative characterization of 
interior voids and surface pockets of proteins [33-35]. Figure 7b shows the binding site pockets in 
the 0-Mad2 state and the intermediate, which are automatically identified by the CASTp analysis. 
The two topological equivalent cavities for the binding site pockets in the 0-Mad2 state and the 
intermediate have different areas and volumes of 183.6 A^/146.0 and 752.6 A^/735.3 A\ 
respectively. This result reveals that the size of the binding site pocket in the intermediate is much 
larger than in 0-Mad2, which suggests that the intermediate is a favorable conformation for Cdc20 
protein binding in the mitotic process. However, in the conformational transition from intermediate II to 
C-Mad2, the distance between the P6/4/5 sheet and the aC helix increases from 15.0 to 19.8 A, which 
increases the locating space for the newly formed P8"/8' sheet and favors the P7/8 translocation to form 
the p8"/8' sheet. 

Figure 6, Mass center distances between the P6/4/5 sheet and aC helix (green), between 
the P7/8 sheet and aC helix (blue), and between the pi strand and aA helix (red) along the 
transition pathway from 0-Mad2 to C-Mad2. The roseate letters included in the grey areas 
denote the different conformations. 
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Figure 7. (a) Moving distance of the P6/4/5 sheet between 0-Mad2 and II; 
and (b) representation of the pockets detected by the CASTp program (cyan spheres: 
Cdc20 binding site pocket). 
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0-Mad2 II 
Total Cavity volume (A' ) 146.0 735.3 

Total Ca\-it>' Area (A- ) 183.6 752.6 (b) 

3.2. Hydrogen Bond and Hydrophobic Structure Analyses along the Transition Pathway from 0-Mad2 
to C-Mad2 

To gain information of hydrogen bonds and the hydrophobic core of Mad2, analyses of hydrogen 
bonds and hydrophobic contacts in the conversion process of 0-Mad2 to C-Mad2 were performed. 
The occurrence of all possible hydrogen bonds for the pi strand and the P7/8 sheet in the conversion 
process are shown in Table 1 by calculating the percentages of times on the corresponding CMD 
simulations. Based on the experimental observations, the hydrophobic core involves the residues 
Leul3 in pi, Phe23 and Phe24 in aA, Phel02 in p5, and Vall93 in p8" [25]. Figure 8 shows the 
hydrophobic contacts at the hydrophobic core. A hydrophobic contact is defined as a distance between 
carbon atoms shorter than 4.5 A. For the conformational transition from 0-Mad2 to II, based on the 
unfolding of the pi strand and the P7/8 sheet in the step of 0-Mad2 to I (as discussed in Section 2.2), 
some hydrogen bonds that initially exist in the 0-Mad2 state disappear gradually with the 
conformational change to I. For example, the hydrogen bond between the N-H group of AsplOl in P5 
and the O atom of Gln9 in pi is maintained with the occupancy of 98% of simulation times in the 
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0-Mad2 state, while the same hydrogen bond is maintained with the occupancy of -0% of simulation 
times in conformation I. Moreover, visual analysis shows that the hydrophobic contact distances of 
3.2, 5.0 and 4.1 A between Leul3 and Phe23, between Leul3 and Phe24, and between Leul3 and 
Phel02 in the 0-Mad2 state increase to 11.5, 7.8 and 12.9 A, respectively, which predicts that 
hydrophobic interaction destruction occurs because Leul3 moves away during pi unfolding 
(see Figure 8). In the formation of the new helix at the iV-terminus of the aA helix in transition 
from I to II, as discussed, some hydrogen bonds are formed with the conformational change to II. 
For example, the hydrogen bond between the N-H group of Serl6 and the O atom of Thrl2 is 
maintained with the occupancy of 98.8% of simulation times in intermediate II, while this hydrogen 
bond is maintained with the occupancy of ~0% of simulation times in conformation I. 



Table 1, The occupancies (%) of hydrogen bonds for the 0-Mad2 state, I, II, III, and 
C-Mad2 state. 



Hydrogen bond 
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II 


Hydrogen bond 
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C-Mad2 


(Ilell)N-H-' 


O(GlnlOl) 


98.2 


0.0 




(Arg99)N-H---0(Alal98) 


0.0 
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• 0(0^9) 
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0.0 




(Arg99)0 •••H-N(Alal98) 
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(GlnlOl)N-H- ■ ■0(Metl96) 
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(Phel51)N-H- 
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0.0 




(GlnlOl)O- ■ •H-N(Metl96) 
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98.8 
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(Asp 1 03)0 • • • H-N(Asn 1 94) 
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(Thrl57)N-H-- 


•0(Hisl91) 


89.8 
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■0(Leul3) 
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■0(Argl4) 
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(Thrl88)N-H •■■0(Hisl91) 
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■0(Glyl5) 
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(Lys200)N-H- • •0(Leu97) 
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(Val20)N-H-- 
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(Ala21)N-H-- 
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0.0 


54.5 









For the conformational transition from II to C-Mad2 in the formation of the new P8"/8' sheet in the 
transition of III to the C-Mad2 state, as discussed in Section 2.2, some new hydrogen bonds are formed 
with the conformational change to the C-Mad2 state. For example, the hydrogen bond of 
NH(Aspl03)-O(Asnl94) in the p8"/8' sheet with the occupancy of 94.2% of simulation times appears 
in the C-Mad2 state. Moreover, visual analysis shows that the hydrophobic contact distances of 13.7, 
10.0, and 12.4 A between Vall93 and Phe23, between Vall93 and Phe24, and between Vall93 and 
Phel02 decreases to 4.2, 4.3, and 4.1 A, respectively, in the C-Mad2 state, which predicts that a new 
hydrophobic interaction is formed with Vail 93 in the newly formed P8" moving into this hydrophobic 
core (see Figure 8). 
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Figure 8. Representation of the hydrophobic interaction in the conformational transition 
pathway from 0-Mad2 to C-Mad2. 




3.3. Alios teric Network Analysis in the Conversion of 0-Mad2 to C-Mad2 

To fiilly understand the allosteric communication in the conversion of the 0-Mad2 state to the 
C-Mad2 state, we explored the relationship between the motion correlations and their conformational 
changes. It can be seen from the motion correlation analysis in Figure 5 that the pi strand correlates to 
the P5-aC loop with large cross-correlate coefficients. Simultaneously, the pS-aC loop correlates to the 
P6/4/5 sheet; then the P6/4/5 sheet flirther correlates to the P7/8 sheet. Such communication transition 
represents the allosteric communication network from the pi strand to the P7/8 sheet through the 
P5-aC loop and P6/4/5 sheet. A coarse-grained picture of correlated segments of 0-Mad2 and C-Mad2 
is shown in Figure 9 with correlation connections. The width of the black-lines connecting correlated 
segments is proportional to the sum of the edges connecting them. The red-lines with arrowheads 
represent the correlation connections. It can be seen that the segment of the pi strand in the 0-Mad2 
state merges into the new segment of pi + aA in the C-Mad2 state, which induces the fluctuation of 
the P5-aC loop due to the fact that the pi strand traverses to the A^- terminus of the aA helix through the 
P5-aC loop. Such allosteric motion produces the correlation connection between the pi strand and the 
P5-aC loop. Because the pS strand and the aC helix are connected to the P6 strand and the P4 strand, 
respectively, the movement of the P6/4/5 sheet to the aC helix relates definitely to the P5-aC loop, as 
expected for the correlations of the P5-aC loop to the P6/4/5 sheet. Furthermore, the motion of the 
P6/4/5 sheet facilitates the translocation of the P7/8 sheet due to the connection of the P6/4/5 sheet to 
the py/S sheet, subsequently causing allosteric communication to the PV/S sheet. 
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Figure 9. Color coded optimal correlated segment network of 0-Mad2 (left) and C-Mad2 
(right). Widths of the connecting black-lines are proportional to the intersegment edges. 
The red-lines with arrowheads and the red dot represent the correlation motion connections 
and the mass center of the P5-aC loop, respectively. 



0-Mad2 C-Mad2 




4. Models and Methods 

4.1. Conventional Molecular Dynamics Simulation 

All CMD simulations for these built models were carried out using the AMBER 9 package [36] and 
froS all atom force field parameters [37-39]. Five Na^ ions were used to neutralize each of the 
0-Mad2 and C-Mad2 models, and an ionic strength of 50 mM was generated by adding 9 Na^ and 
9 Cr ions for these two models [14]. Similar counterion processes were applied to the other models. 
Each of the systems was explicitly solvated by using the TIP3P water potential inside an orthorhombic 
box of water molecules with a minimum solute-wall distance of 10 A. The protocol for all CMD 
simulations is described as follows: (1) the systems were energetically minimized to remove 
unfavorable contacts. Four cycles of minimizations were performed with 5000 steps of each 
minimization and harmonic restraints on the Mad2 protein from 100, 75, 50 to 25 kcal/(mol-A^), which 
means that the restraints were relaxed stepwise by 25 kcal/(mol-A^) per cycle. The fifth cycle consists 
of 10,000 steps of unrestrained minimization before the heating process. The cutoff distance used for 
the non-bonded interactions was 10 A. The SHAKE algorithm [40] was used to restrain the bonds 
containing hydrogen atoms; (2) Each energy-minimized structure was heated over 120 ps from 0 to 
300 K (with a temperature coupling of 0.2 ps), while the positions of Mad2 protein were restrained 
with a small value of 25 kcal/(mol-A^). The constant volume was maintained during the processes; 
(3) The unrestrained equilibration of 200 ps with constant pressure and temperature conditions was 
carried out for each system. The temperature and pressure were allowed to fluctuate around 300 K and 
1 bar, respectively, with the corresponding coupling of 0.2 ps. For each simulation, an integration step 
of 2 fs was used; and (4) Finally, conventional molecular dynamics (CMD) runs of 50 ns for the 
C-Mad2 and the 0-Mad2 models were carried out, respectively, by following the same protocol. 
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4.2. Targeted Molecular Dynamics Simulation 

Targeted MD simulation [41,42] is a method to observe large-scale conformational transitions 
between two known end-point conformations of a molecule. A restraint energy term was added to the 
energy function proportional to the square of the difference, which may be characterized as the 
mass-weighted root-mean-square deviation (RMSD) of the current structure to the target structure in 
terms of atomic positions. The fiinctional form of the restraint energy can be written as: 

^TMD = ^kN [RMSD it) - RMSDoit)f (1) 

where, k is the harmonic force constant per atom, is the number of the restrained atoms, RMSD (t) is 
the root-mean-square deviation of the simulated structure at time t relative to the target structure, and 
RMSDq {t) is the prescribed target RMSD value at time t that decreases to zero linearly with time to 
drive the system from an initial structure to the target structure. To obtain the appropriate small force 
constant k, four independent, short (4 ns) TMD simulations were performed using different force 
constants of 0.2, 0.5, 0.8 and 1.0 kcal/(mol-A^), respectively. As shown in Figure S3 for the conversion 
of 0-Mad2 to C-Mad2, three out of four simulations reached a RMSD value of <2.0 A extracting from 
the backbone atoms of the Mad2 protein. Therefore, k = 0.5 kcal/(mol-A^) was chosen as the lowest 
harmonic force constant to apply onto all the backbone atoms of the Mad2 protein to bias the 
trajectories toward the target structure. All the other atoms in the Mad2 protein were assumed to move 
freely and to rearrange at their convenience to let the protein reach its equilibrium under this 
constraint. The weighted RMSDs along the TMD trajectories from four, 4 ns independent simulations 
showed that each process could follow the same conformational transition pathway, which indicated 
that the external forces and initial velocities did not affect the transition pathway. Therefore, the final 
TMD simulations were performed with the 10 ns transition TMD, followed by 20 ns of equilibration 
with an integration step of 2 ps. To test the mutual conversion property between the 0-Mad2 and 
C-Mad2 states, the TMD simulations were performed in an 0-Mad2-«-^C-Mad2 bidirectional fashion. 

4.3. Correlation of Atomic Motions 

The dynamic feature of a protein and the extent of correlation of motions in the different regions of 
a protein were assessed via the calculation of cross-correlation coefficients, C{ij), as follows: 

C(i,j) = (An X Ar;/[ (Ari2)(Ar/)]i/2 (2) 

In the equation, Ar, and Ar, are the displacement vectors for Ca atoms of residues i and j, 
respectively, and the angle brackets denote the ensemble averages [36,43]. In the present study, the 
correlation coefficients were averaged over the regions of protein, and the resulting cross-correlation 
coefficients are presented in the form of a two-dimensional graph. These structure analyses in the 
present work were calculated by using the PTRAJ module of the AMBER 9 program. 

4.4. Surface Analyses 

Shape descriptors representing protein structure, such as depth, surface curvature, extreme 
elevation, surface area, and volume, have been used extensively to identify, study, and compare 
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protein-ligand interactions, protein-protein interactions, and the respective binding sites. 
We used the Computed Atlas of Surface Topography of proteins (CASTp) program 
(http://sts.bioengr.uic.edu/castp/) [44] to determine the binding pocket. The CASTp program uses 
weighted Delaunay triangulation and the alpha complex for shape measurements. It provides 
identification and measurements of surface accessible pockets as well as interior inaccessible cavities 
for proteins and other molecules. It analytically measures the area (in A^) and volume (in A^) of each 
pocket and cavity, both in solvent accessible surface (SA, Richards' surface) and molecular surface 
(MS, Connolly's surface) [44]. In our work, the probe radius was set to the default value (1.4 A). 
The detected pockets from using these algorithms are ranked with their volumes and areas. 
For ligand molecules, internal protein cavities appear to be a favored binding site, and the cavity 
volume may play an important role in the strength of the guest molecule-host cavity interaction. 
In CASTp analysis, we chose the Cdc20-binding pockets for our analysis. 

5, Conclusions 

Conventional molecular dynamics simulations and targeted molecular dynamics simulations were 
performed to address the stable structure characteristics for the open and closed Mad2 states and the 
dynamic transition mechanism of the two states, respectively. The results of this study show that the 
predominant structural differences between the stable open and closed Mad2 states were the 
transformations of the P7/8 sheet and the pi strand with the rearrangement of P7-P8-P6-p4-P5-pi to 
P6-P4-P5-P8"-P8'. The results for the transition mechanism investigation demonstrated that the 
structural transition from the open Mad2 state to the closed one could occur via an intermediate. 
The structural change characteristics along the transition pathway indicate that the pi strand and the 
P7/8 sheet unfold first, then move, and finally form an additional helix and a new P8"/8' sheet via the 
intermediate. Particularly, the pi strand traversing through the P5-aC loop is a prerequisite for the P7/8 
translocation. The conformation analysis for the intermediate indicates that the movement of the P6/4/5 
sheet close to the aC helix causes the exposure of the P6 strand and an increase in size of the active 
pocket, and favors the binding of the Mad2 protein to Cdc20. This observation reveals the critical 
function of interconversion of the two states in APC/C inhibition. Correlation analysis reveals that the 
allosteric network occurs from the pi strand to the P7/8 sheet via the P5-aC loop and the P6/4/5 sheet 
in the conversion of the open Mad2 state to the closed Mad2 state. The present investigation provides 
useful insights into understanding the dynamics of the transition mechanism of the active closed Mad2 
state to the open Mad2 state. 
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